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For quantum systems with two dimensional configuration space we construct a physical radial 
momentum observable. Rescaling the radius we find the dilatonic degrees of freedom form a Weyl 
àlgebra. With this we construct the radial Wigner quasi-probability distribution function. 
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The Wigner quasi-probability distribution function is a 
familiar tool to many working in quantum and atom òp- 
tics jp . It is primarily used in the classical-quantum cor- 
respondence where the appearance of positive and neg- 
ative regions of the Wigner function gives easily under- 
stood information concerning the probability concentra- 
tions and quantum intcrfcrcnces present within the quan- 
tum state g. Typically one describes the Wigner func- 
tion on a phase space which is labelled by cartesian coor- 
dinates of position and momentum. For physical systems 
which admit a two dimensional cylindrical symmetry, eg. 
trapped ultra-cold ions, bose condensates, etc, clearly a 
polar description of the Wigner function would be more 
natural. However, no such description has appeared in 
the literature. In this letter we show, in three stages, 
how a Wigner function for the radial observables can be 
constructed. This radial Wigner function could be recon- 
structed from experimental data much as recent experi- 
ments have reconstructed the cartesian Wigner function 
for the one dimensional motion of a trapped ion || . The 
angular parts of the complete four dimensional Wigner 
function are complicated by the imposition of singleval- 
uedness of the wavefunction under a rotation of 2tt which 
cause the conjugate angular momentum to become dis- 
crete. We leave the angular part for a later work. 

The stages towards the construction of a radial- Wigner 
function proceed as follows: (1) a proper Wigner func- 
tion possesses marginals which are true probability dis- 
tributions for the observables whose eigenvalues label the 
Wigner function and thus the phase space axes. For 
a single dcgrcc of freedom a mere transformation of 
the cartesian position and momentum into polar form 
does not yield a proper Wigner function for the po- 
lar observables This is also true for higher dimen- 
sional phase space representations. Central to the prob- 
lem is the correct specification of the radial momentum 
operator. By noting the symmetry action of the mo- 
mentum on the half-infinite radius observable we con- 
struct a physical "conjugate" momentum P r ; (2) essen- 
tial to the construction of the Wigner function is the 
existence of "point" operators A(Ai,A2), which obey, 
Tr[i(A 1 ,A 2 )it(A' 1 ,A^)] = «S(Ai - X[)S(X 2 - A 2 ) 10. 
We find that the radial A(r, P r ), is not completely expo- 
nential in the radial position and momentum operators; 
(3) guided by the form of A(r,P r ) we make the coor- 
dinate transformation v = lnf. We find, [w,-P r ] = ift, 



the Weyl àlgebra. After re-scaling the eigen-basis kets 
of v to recover the Standard resolution of unity we can 
easily construct the radial Wigner function, W(v,P r ). 
This radial (or dilation) Wigner function gives the proper 
marginal probability distributions for v and P r . We also 
note the existence of a dilaton ground state |0) düations, 
and calculate the wavefunction of this ground state in the 
f basis. We finally calculate the radial Wigner function 
for the lowest Schwinger states |í,0), (these are simul- 
taneous eigenstates of energy and angular momentum), 
and briefly outline how, given a quantum state in a two 
dimensional harmònic Fock representation, one can con- 
struct the radial reduced density matrix, and from there 
the radial- Wigner function. 

Dirac in his textbook on quantum mechanics ^| intro- 
duced the following momentum, conjugate to the radial 
coordinatc 



P D = \ (xP x + yP y - ih/2 



(1) 



The factor of two difference in the U term arrises because 
we are working in two dimensions instead of the three in 
Dirac's case. Using the cartesian commutation relations 
between x and P x one can easily confirm that 



(2) 



Using the Baker-Campbell-Hausdorff (CBH) expansion 
and (121) we can show 



exp(içP u /h) f exp(-içP u /%) = f + ç 



(3) 



Since f has a half-infinite spectrum, the operator P D 
cannot be self-adjoint. Although this has been noted by 
many others no self-adjoint radial momentum has 
been proposed that the author can find. Before going on 
to construct a self-adjoint radial momentum we first cal- 
culate the representation of P D in the f eigen-ket basis 
(i.e. f\r) = r\r)). From the commutation relation (||) we 
must have (r\P D \ip) ~ —ih(d r + f(r)){r\ip). From (Q), 
and expressing the cartesian partial derivatives in polar 
coordinates we find 



(r\P D \Íj) =-in[d r + ± ) ;r|,-> 



(4) 



To obtain a self-adjoint momentum operator we must find 
an operator which respeets the half-infinite spectrum of 
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f. This can be achieved if we have [f , P r ] — iUf, the 
Sack dilation àlgebra ||. Using the CBH expansion we 
can easily show 



D(X,fj,) = exp(ifiP r /2%) r lX exp(ifiP r /2h) 



(10) 



exp(içP r /?i) f exp(— içP r /Ti) = fe} 
To explicitly construct a P r we form 



(5) 



(0) 



where the additional h/2 makes P r hermitian. To obtain 
the representation of P r in the \r) basis we simply use 
(f|) above to find 

(r\P r \ip) = -m(rd r + l)(r\ij) . 

From this we can calculate the transition function (r\P r ) 
where \P r ) is the eigen-ket of the operator P r to be 

l r iP r /h 



(r\P r ) = -== . (7) 



/2wh r 

This is normalised with the measure rdr and gives 

/ rdr (r\P r )(P' r \r) = S (P r - P' r ) . 
o 

The action of f on P r is not a simple scaling, instead 

e -içf/hpr e içf/h _ pr _|_ - 

To form a displacement, or point, operator we might be 
tempted to exponentiate a linear combination of the po- 
sition and momentum operators in analogy with the har- 
mònic oscillator. However this construction does not re- 
sult in a unique adjoint action of the displacement oper- 
ator on f and P r . From [|10| we see 



T>(a, m) = exp 



la 



exp 
exp 



h 

ia 

¥ 

imf 



P r + rr 
P r ^j exp 
(e a 1) 



(8) 



imr 
exp 



(1 - e -) 



-P' 



Thus the choice of ordering alters the action of T>(a,m) 
on the operators f and P r . This is to be expected since 
the àlgebra is not symmetric in f and P r . One must ex- 
plicitly specify the ordering in this displacement operator 
as is done in most treatments of displacement operators 
for spin systems. 

An essential property for a displacement operator is 

10. 



Tr[D(X, n)D\\', //)] w S(X - X')S(fx - //) 



(9) 



No operator in the form of products of exponentials of f 
and P r will satisfy (^). We instead consider the operator 

lE3, 



where À G [0, +oo). Taking the trace over |r) with mea- 
sure íp-, from r — to r = +oo gives 
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Tr[D(X 7l· L)Di(X', l· i')} = -^S(X-X')ó( t ,- t /) . (11) 



One can show D'(X,fi)fD(X ) ^,) — fe~ M . However to 
evaluate D'(X, fj,)P r D(X, /i) one must calculate f a p r f~ a . 
This is done by expressing P r in terms of P D through 
(^) , using the simpler commutation relations for P D , the 
rule [f(f),P D ] — ihd f (f) / df , and finally transforming 
back to P r . This surprisingly gives D' (A, [i)P r D(X, fx) = 
P r + X. 

Thus, by using the new "displacement" operator (|k|), 
one can (except for the exponential factor in (|ïï|)), almost 
completely regain the properties of the Harmònic oscil- 
lator displacement operator. Although one could now 
proceed to construct a Wigner function with this dis- 
placement operator, the replacement of expíiAf) with f lA 
in (|Ï0|), and the exponential prefactor in ( |lï| ) point to- 
wards a clearer understanding of the situation. Guided 
by these signals we now make the transformation to a 
new coordinate operator, 



lnf 



(12) 



Using the defmition of P r in terms of P D (||), we can 
show 



[v,P r ] = [lnf,P r ] = [\nf,fP D ] = f[lnf,P D ] =f—=íh 

f 



(13) 



With this we have re-gained the àlgebra of the original 
cartesian position and momentum observables. Since the 
spectrum for both v and P r ranges from (— oo, +00) the 
translation generating Weyl àlgebra in ( |l3| ) properly re- 
spects the spectrum of both operators with the result 
that both v and P r are self-adjoint and represent physical 
observables. We can construct creation and annihilation 
operators in the same manner as the harmònic oscillator, 
à = (v + iP r ) /-v/2, and fit = ({) _ iP r )/y/2. In particular, 
there exists a proper vacuum state |0) düations which is 
annihilatcd by à. 

Not everything is straightforward however. The reso- 
lution of unity in the |r) basis becomes, in the eigenbasis 
of v, 



+00 />oo 

dx \x)(x\ — / rdr \r)(r\ 
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= / e Zv dv \v)(v\ 
With this resolution of unity we must have 



(14) 
(15) 
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TL\v') = \v') = 



2vr 



e 2v dv \v)(v\v') 



(16) 



W(Ç,a) 



d 2 a 



C(a,s) 



(19) 



and thus (v \v') = -^5{v — v') = 25{e 2v — e 2v ). This ex- where C(a, s) is the s-ordered displacement operator 



ponential measure complicatcs mattcrs and so to simplify 
we rescale the basis kets to be 

\v)' = e^) r (17) 

where \e v ) r is a eigenket of f with eigenvalue e v . The 
resolution of the identity in the new l'ü)* basis is, 27rl = 

The inner product between \v)' basis 
kets is now simply *(v\v') m = S(v — v'). This rescaling 
will primarily appear in the calculation of '(v\ip) when 
we have (r\ip). To be certain that everything is identi- 
cal to the Harmònic oscillator case we finally evaluate 
* {v\P r \ip) . Since \v)* = r\r) where r = e v , we have 



'(v\P r \ip) = r{r\P r \iP)\ 
-ih(rd r + 1) 



r—e 
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r[-ih(rd r + l)] (r#)| r= 



= -iKrd r ' (v\ip)\ r=e v = -ihdv' (v\ip) . 

We see that the momentum operator P r , in the basis \v)* 
is the familiar —ihdy, the same as the cartesian case. 

As in the Harmònic oscillator case we define the 
Wigner and s-ordered quasi-distribution function to be 



C(a, s) = Tr pD(a)e 



M 2 /2 



(20) 



and D{a) = exp(aà j( — a*a), where and à , the dila- 
tonic creation and annihilation operators, are defined as 
above. For s = 0, £ = 7 + iS, the Wigner function be- 
comes 



win, 8) 



1 

27 



de e 



-ieS 



(21) 



We now calculate the Wigner function for the quantum 
states pi = \l, 0){l, 0|, where \l,m) is the Schwinger angu- 
momentum state given by 



\l, m) 



y/{l + m)\{l-m)\ 



-,A 



\{l+m) ít(i-m) 



0,0) , (22) 



where À + = (à x — ià y )/-\/2, À- = (à x +ià y )/\/2 Ei- 
ther by solving the generating diffcrcntial cquation (as is 
done in flÏ3| ) or by solving the harmònic oscillator radial 
equation one finds (r, <t>\l,m) — Ri, m (r)&m(<P), where 



Rl 



i(r) = ^\Zfrw (/3r)2|ml e ~" 2r2/2 l hh (/3V) ( - 1)Hm| ' 



e m (0) 



2in%4> 



2tï 



(23) 
(24) 



L" is the associated Laguerre polynomial and (3 — l/yh. Integrating out the 4> dependence, setting h = 1 and using 
( |l7| ) we obtain 



m (v\l,0) = v / 2e s e- e2B/2 L?(e 2s )(-l) í 

Inserting (^5|) into ( pï| ) we finally get 

2e 2 7 r+oc 



n 27 r-t-oo 

Wi(j,S) = — de e- 2íeí exp [-e 27 cosh2e] L i (e 2(7+£) )L i (e 2( ^- £) ) 

^ J — OC 



(25) 



(26) 



The radial Wigner function, ((2^), is numerically in- 
tegrated for the first four \l, 0) states and plotted in 
Figs. l(a)-l(d). The presence of negative regions in 
these Wigner functions is not surprising as the pure 
states \l, 0) are reminiscent of Fock states. One can con- 
struct the dilatonic coherent states \a) = D (a)\0) düatons- 
More interesting is (r\0) düatons ■ From the definition 
à\0)düatons = 0, an d \/2à = v + iP r , we can evaluate 
' (v\à\0) düatons = to get * (v 1 0) düatons = N exp(-v 2 /2), 



where iV 2 = y/ïr. Again using (Jïï]) we obtain 

1 l r -í lnr 

(r\0) düatons = — = exp(--ln?Tnr) = -=r- , (27) 

r^ir 2 ry7r 

which is normalised to unity, f Q °°rdr | (r |0) düatons | 2 = 
1. 

In this letter we examined the problem of constructing 
a Wigner function for the two dimensional radial sub- 
space of a quantum system possessing two continuous 
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dcgrccs of freedom (or a four dimensional Wigner func- 
tion). Since thc radial subspace is labelled by an op- 
erator with an half-infinite spectrum, previous attempts 
to define a physical conjugate momentum have failed. 
By choosing a momentum operator which respected the 
spectrum of the radial coordinate we found we could con- 
struct a physically meaningful radial conjugate momen- 
tum. With a logarithmic rclabelling of the radial co- 
ordinate the new coordinate and momentum satisfied a 
Weyl àlgebra and we were thus able to carry over all 
the techniques associated with the harmònic oscillator to 
construct a Wigner function. One can identify a ground 
state for these observables which is destroyed by a dila- 
tonic annihilation operator. We finally examined the 
radial- Wigner function for particular quantum states and 
found the radial wavefunction for the dilatonic vacuum 
state. To examine the radial- Wigner function for more 
general quantum states is difficult. If one has the den- 
sity matrix in a x — y Fock state basis, one must first 
transform p into the Schwinger basis of angular momen- 
tum kets, p = ^2 Ci m: i' m i\l, m){V ', m'\. One can then find 
(r,(f>\p\r' ,<fi'). Finally one must trace over <f>. However, 
even though we have traced out over <j> one will still end 
up with a sum over m. This is so because, in contrast 
with the cartesian Fock state decomposition, the radial 
wavefunction is labelled by both quantum numbers l and 
to while the angular part is only labelled by m. We 
also note that the Wigner function here gives the cor- 
rect marginals for the operators P r and v. Once these 
marginals and their corresponding integration measures 
are obtained one can rescale v to r. One cannot rescale 
the axis in Fig. 1 to find the Wigner function for P r 
and f without altering the integration measure and thus, 
the function form of the Wigner function. Finally, al- 
though an observable, the radial operator v — lnf, may 
not be easy to measure directly. However, its should be 
possible, through Standard reconstruction techniques, to 
numerically approximate the radial Wigner function from 
experimental data. 
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FIG. 1. (a)-(d) Radial Wigner Functions for the states 
|í,0>, í = 0,..,3. 
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